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C^! Abstract 

D I The tetramixing of pseudoscalar mesons Tr-rj-rj' -rjc and vector mesons uj-p-c/^-J/ip are studied in 



^ 



the light-cone constituent quark model, and such mixing of four mesons provides a natural source 



^ , for the intrinsic charm cc components of light mesons. By mixing with the light mesons, the 

OO , 

QQ ' charmonium states J/t/j and rjc could decay into light mesons more naturally, without introducing 

o 

lO ' gluons or a virtual photon as intermediate states. Thus, the introduction of light quark components 

f— ^ ■ into J/t/j is helpful to reproduce the new experimental data of J/tp decays. The mixing matrices and 

the Q^ behaviors of the transition form factors are also calculated and compared with experimental 

data. 

X 

C^ ■ PACS numbers: 12.39.Ki, 13.20.Gd, 13.40.Gp, 14.40.Be 



* Electronic address: mabq@pku.edu.cn 



The mixing of mesons has been widely investigated since the 1960s, when the concept 
of a mixing state of the p-oj mesons was proposed [l| by considering that the electromag- 
netic interaction does not conserve isospin. Later, the u-ip mixing and rj-rj' mixing were 



introducec 



formula 



2| to explain the deviation of the meson mass from the Gell-Mann-Okubo mass 
^. It was also pointed out that the difference between the u and d quark masses 



7| and p-u- 



aa 



were 



introduces the Ti-rj mixing [5|. Then, the trimixing of 7r-r]-r]' [&, 
proposed, and their effects were studied in different methods. On the other hand, the cc 
contribution to the r] and r]' mesons was considered |9|, and the trimixing r]-r]'-r]c was stud- 
ied 10|. As a further extension in this paper, we try to combine the above two types of 



trimixing by considering the tetramixing of pseudoscalar mesons TT-r]-ri'-r]c. The mixing of 
gluon component gg and r]-r]' were also studied |lll . Il2| . As the mixing of r] and rj' is still 
not completely clear right now, we think that the charm and gluon components may both 
be possible to mix with these mesons, and it is worthwhile to study both of them carefully. 
The recent CLEO experiment 131] of the charmonium decays J/ip — )■ 77r, 777, and 777' 
also motivates us to extend our tetramixing to the vector mesons uj-p-(j)-J /il). According to 
the pure valence cc structure of charmonia in the naive quark model, these decay modes of 
charmonium ipijiS) must happen via the annihilation of the heavy quark constituents into 



gluons or a virtual photon [l3|,[lj], because of the Okubo-Zweig-Iizuka rule, which postulates 
a suppression of transitions between hadrons without valence quarks in common [ij]. More- 
over, the mechanisms of these decays are not completely clear yet, and there are various ways 
to describe them, such as ip{nS) — )■ '-ygg — )■ 7P, ip{nS) — )■ ggg — )■ gg7, ip^nS) — )■ 7* — )■ qqj, 
and so on ISj. However, with the model of uj-p-(j)-J /ip mixing, the above-mentioned decays 
of J /ip could occur more naturally through the direct transition from its light quark com- 
ponents to light mesons such as vr, r], or rj' without introducing intermediate gluons or a 
virtual photon, and the cc components of these light pseudoscalar mesons also allow J /ip to 
decay to them. The mixing of u-p-cp-J /ip is thus helpful to reproduce the new experimental 
data of J /ip decays. 

For the light vector mesons w, p, and 0, the existence of cc states in them may be 
interpreted as a support to the theory of intrinsic charm [l5| in these mesons. Different from 
the extrinsic quarks, which are generated on a short time scale in a reaction process with 



large momentum transfers, the intrinsic quarks are intrinsic nonperturbatively to t 



le hadron 



wave function and exist over a time scale independent of any probe momentum 15|, |l6| . The 



postulation of intrinsic cc components in p and vr offers a possible solution of the "pvr puzzle" 
by allowing direct transitions between J /ipi^p') andp(7r) through the rearrangement of the 
valence and the intrinsic cc components of p(7r) pM\- Now the tetramixing of u-p-cjy-J /if} 
introduces the intrinsic cc components into all three light vector mesons w, p, and 0, and 
J /ip can decay to them in a similar way, without annihilation of the quark constituents. 
This applies to the pseudoscalar mesons vr, 77, and r^', too, as they mix with the charmonium 
Tjc in our model. The intrinsic cc component in rj' was also studied in Refs. 17|, [iSj, in which 
the intrinsic charm content of the rj' meson fZ was evaluated, and we shall compare our 
result of /5 with previous results in Refs. 17|, ll8|] and other works at the end of this paper. 



We adopt the light-cone constituent quark model 19|-|21| to study the mixing of mesons. 
The light-cone constituent quark model is a convenient and effective model to treat the 



nonperturbative aspect of QCD, and the mixing of mesons in this model has been studied 22, 
23|. 

The mixing of pseudoscalar mesons and vector mesons could be described by two SO (4) 
rotation matrices My and Mg, respectively: 
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in which the unmixed states are 



^i = -jniuu + dd)ip^j, pi = -js{uu - dd)ippj, (pi = -ssip^^, J/ipi = ccipj/^^, 



V2 



(2) 



TTi 



75 (mm - dd)ip^j, 7]q = ^{uu + dd)ipr,,, Vs = ssipns, Vco = cc^vcoi 

where y? is the momentum space wave function of the corresponding meson. 

Since the rotation group SO {4) = 5*0(3) ® 5*0(3), the 5*0(4) mixing matrix M can be 
written as 

M = R+R_, (3) 



where the matrices i?+ and -R_ are generated by the 5*0(3) generators, and each of them 



could be parameterized by three independent rotation angles as 



R, 



( 



n-, 



02, O3) 



cos ^ 
^sinf 

a 2 

^sinf 

a 2 

\ ^ sin f 

\ a 2 



-^sin^ -^sin# -^sin^\ 



cos I 
^sinf 

a 2 

-^sinf 

a 2 



sm 

-^sinf 

a 2 



COS ■ 



^sinf 

a 2 



■ sm ■ 

ut z, 

fesinf 

a 2 



— ^ sm 77 

a 2 



cosf / 



(4) 



where 



R. 



Of, , Op. 



^4, '^51 



COS ■ 



.64 sin £ 

/3 ^"^ 2 



^sin^ 

13 ^^" 2 



-sm| 



|sinf \ 



cos I -^sm| ^sm| 



2 /3 



COS ^ 



|sinf -|sinf |sinf 



-^sin^ 
cos 4 



a 



Ol + Ol + Ol (5 = ^01 + 01 + 01 
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and thus the mixing matrix M is parameterized as six independent rotation angles (^1, ^2, 
... ,6*6). Our detailed procedure of obtaining the matrix form of -R+ and -R_ [Eqs.(4) and 
(5)] is given in Appendix 1X1 

When referring to the mixing of specific types of mesons, M stands for My or M^, and 
the parameters change to {0\, O2, ...) or {Of, 6^1, ...) correspondingly. 

During the numerical calculation, we also used a more compact form of M with eight 
real parameters under constraints, and the detailed procedure is given in Appendix [Bl 



The decay constants and transition form factors also mix as 
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The above decay constants and transition form factors are defined as 23|, |24| 



{0\j,\V{p,S,)) = Mvfve,{S,), 

{lip - q)\J'\P{p, A)) = te'Fp^,,,{Q')e^^^^p,e,{p - q, A)g., 
{P{p')\J^\V{p, A)) = teFv^P,{Q')e^'''''e,{p, X)p'p,. 



(8) 



(9) 
(10) 
(11) 



To calculate them, we use the light-cone quark model with the Fock state expansions of the 
unmixed mesons [the right-hand side of Eq. ([1])]: 



1^) = ^\(i(i)i^qq + ^\(i(i9)i^qqg + 



(12) 



and to simplify the problem, we adopt the lowest order of the above expansions, which takes 
only the quark- ant iquark valence states of the unmixed mesons into consideration. 
The wave function of an unmixed meson in the light-cone formalism is 

dxd^k I 



19 



25| 
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(13) 



where y? is the momentum space wave function, described by the Brodsky-Huang-Lepage 
prescription |l9l. |25|: 



(p{x,kj_) = ipBRL{x,k±) = Aexp 
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(14) 



[A and /3 are the parameters of the meson, and mi and 7112 are masses of the constituent 
quarks), and x^|(x, k_L, Ai, A2) is the hght-cone spin wave function, which is related to the 
instant-form spin wave function by the Melosh-Wigner rotation 261428 1 



xUt) = w,[{kl + m,)x]{F) - A;fxt(F)]; x\{T) = w,[{kl + m,)x\{F) + fcf xl(i^)], (15) 



where Wj = 1/ ^2kl{kS' + m^), k 



R,L 



k^ ± k^, k+ = k^ + e 



xM, and M 



k_L +m.1)/x + {'k± +m.2)/{l — x). The Melosh-Wigner rotation is an important in- 



gredient of th e lig 
spin puzzle" [28, 



it-cone quark model and plays an essential role in explaining the "proton 



29| . The detailed formulas for calculating the decay constants and tran- 



sition form factors of mesons were listed in Ref. 23|, and the examples of applying them 
to set meson parameters and to calculate the decay constants and transition form factors 
numerically can be found in Ref. 30 1. 

The values of the meson parameters A, /3, mi, and m,2 and the parameters of the mixing 
matrices (a„, b^, ...) and (a^, b^, ...) can be chosen by fitting the light-cone constituent 
quark model results of the meson decay constants and transition form factors (at Q^ = 0) 
to experimental data. The Q"^ -^ 00 limiting behavior of Q'^Fp^^^* is also considered as a 



constraint to set the parameters 21 



lim Q'Fp^^^,iQ^)=2cpfp 
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(1, 5/3, a/2/3). With the similar method as Ref. jl7|, we 
^^Q — 4a/2/3. All these requirements are taken as constraints to determine the 



where cp 
also obtain c 

parameters of mesons and parameters of the mixing matrices. During our calculation, we 
first use the decay constants and the radii of tt"*" and K~^ as the constraints to locate the 
values of A^, /3^, m„, and m^, assuming that the wave function parameters of vr^ are the 
same as those of tt/, and the constituent quark mass m,^ ^ m,u (their difference could be 
ignored compared with rric) [7j^. The other parameters are then determined under the left 
constraints. 



Our numerical calculation gives the mixing matrices of vector and pseudoscalar mesons: 

/ 0.9886 -0.0122 -0.1429 0.0076 \ 

0.0299 0.9910 0.1221 -0.0025 

0.1400 -0.1250 0.9808 0.0258 
\ -0.0111 0.0058 -0.0239 0.9986 / 



M„ 



(17) 



M. 



(18) 



/ 0.9895 0.0552 -0.1119 0.0342 \ 

-0.1082 0.8175 -0.5614 -0.0259 

0.0590 0.5696 0.8160 0.0452 

\ -0.0395 -0.0065 -0.0478 0.9960 / 

We see that some of the entries of the mixing matrices are small; for example, one of the 
entries in the first row of My is 0.0076. But this nonzero entry means a charm component 
in the u meson, which allows r/c to decay to u directly in our model. Other entries of the 
mixing matrices have the same meaning, and it is such entries that are helpful to reproduce 
the experimental decay data of J/ip and other meson decays. 

The results of fitting light-cone constituent quark model results to experimental data 
are shown in Table [B The fourth column contains the results of tetra-mixing model ir-r]- 
rj' and p-u-(j), while J/tp and rj^ do not mix with other meson states, with the values of 
their parameters (MeV) set as Aj/^ = 31.1660, /3j/^ = 0.9777, A^^ = 125.7935, and f3rj^ = 
0.7524 to fit the experimental data. The most apparent differences between tetramixing and 
trimixing results are in the last four rows, which show that the trimixing formalism do not 
explain the nonzero decay width of J/ip — > vr, t], and t]', while the tetramixing formalism 
can well reproduce these experimental decay data. The parameters of the mesons and the 
mixing matrices determined during the fitting process are listed in Tables [III and Table lllli 

The Q^ behaviors of the form factors of vr, rj, and t]' are shown in Figs. 1-3, and we 
see that they are generally in agreement with the experimental data. The Q^ behavior of 
the form factor of rjc is shown in Fig. 4, and, by comparing with theoretical data from 
another model, the calculated curve fits well in most of the lower Q^ region. We can also 
obtain the Q^ behavior of the transition form factors in the timelike region, either by making 
the substitution qj_ — )■ iqj_ [33] or by parameterizing the transition form factors as explicit 
functions of g^ in the spacelike region and then extending them through analytic continuum 



TABLE I: Experimental data and the light-cone constituent quark model fitting results of the 



meson decay constants and transition form factors. The experimental data (unmarked) are from 
Ref. [3.4], and the experimental data (marked with daggers) are from Ref. [13']. The data in the 
fourth column (unmarked) are from Ref. [7|, and the data (marked with stars) are calculated 
assuming that J /tp and 77c do not mix. 



Decay Constants 


Experimental 


Theoretical Fitting 


Theoretical Fitting 


or Form Factors 


Data (GeV) 


of tetramixing (GeV) 


of trimixing (GeV) 


Fn^-yy' (0) 
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2.4639 
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0.4494 ± 0.0197 
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Fp^TT-y* (0) 


0.8237 ± 0.0549 
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Fp-^n'y* (0) 


1.5687 ± 0.0525 
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-^»7'-^P7*(0) 


1.3175 ± 0.0327 
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F^^TTJ* (0) 


0.1331 ± 0.0032 


0.1301 


0.132 


F^^^^f* (0) 


-0.6937 ± 0.0071 
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Fff.^n'j* (0) 


0.7153 ± 0.0125 
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Fj/^p^n.r (0) 
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TABLE II: The meson parameters A and /3 (GeV), and the masses (GeV) of constituent quarks 
determined from the fitting process. 

Au}j Apj Afj,j ^J/ipi ^TTi A^^ Ari^ ^r]cO ^^u{d) "T-s 

41.4712 38.1430 63.1638 31.1724 47.3635 38.7860 95.4496 125.8099 0.198 0.556 

l^uji l^pj f^^j I^J/^j f3nj f3r,q f^ris Pr^^o "^c 

0.4319 0.4318 0.4757 0.9781 0.4112 0.4887 0.4887 0.7373 1.270 
TABLE III: Parameters of the mixing matrices My and Mg determined from the fitting process. 
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to the timelike region 3J]. The results are shown in Figs. 5-10, among which Figs. 5 and 
Fig. 6 are compared with the experimental data, while Figs. 7-10 could be considered as our 
predictions of the Q^ behaviors of J /ip transition form factors. 

We can further use our results to learn the properties of the intrinsic cc component in 
the light pseudoscalar mesons. With the Fp_>^^. (0) (where P = tt, rj, rj', t]c) in Table [land 
the mixing matrix Ms, we obtain the transition form factors of unmixed mesons Fp^^^^*{0) 
(where Pj = ttj, rjg, rjs, rjco)- Taking them into Eq. (TT6|) . we have the values of /^, fg, fs, 
and U 



U = 0.0984 GeV, 

(19) 



fg = 0.0976 GeV, 



/, = 0.1298 GeV, 
/e = 0.4874 GeV. 
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FIG. 1: The Q^ behavior of the form factor Q"^ Fj^^^^* {Q"^) compared with experimental data 
431. 



m 



/ f-K fg fs fc \ 

£n £q £s £c 

J rj J r] J rj J r] 

f-rr £q £s £c 

J Tj' J r]' J Tj' J Tj' 

\ r f? f* f^ / 



M. 



-0.0145 0.0168 \ 



(20) 



f U \ 

A 
/, 

V fj 

( 0.0974 0.0054 

-0.0106 0.0798 -0.0729 -0.0127 

0.0058 0.0556 0.1059 0.0219 

\ -0.0039 -0.0006 -0.0062 0.4855 / 

We see that /^, = 0.0219 GeV = 21.9 MeV. It is compared with previous results in TablejIVl 
/5 could be considered as the reflection of the intrinsic charm content of the v( meson jlTJ, 
and we see from Table IIVI that our result of /5 is in the similar region with most of the 
previous results. 

In summary, we use the light-cone constituent quark model to study the tetramixing 
of pseudoscalar mesons vr-r^-r^'-r^c and vector mesons uj-p-(j)-J /ip. The parameters of mixing 
matrices and meson parameters are determined by fitting our theoretical model results of 
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FIG. 2: The Q^ behavior of the form factor Q'^F^^^^^* (Q^) compared with experimental data 42, 
431. 



TABLE IV: The value of f^, (MeV) in different models. 

Our Feldmann Halperin and Cheng and Cao, Cao, Huang Yuan and 
model and Kroll [35] Zhitnitsky [36] Tseng [37] and Ma [17] Chao [38] 
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the meson decay constants and transition form factors (at Q^ = 0) to the experimental 
data. We also calculate the Q^ behaviors of the meson transition form factors, and these 
results are generally in agreement with the experimental data or results from other models. 
Our results of the Q^ behaviors of transition from factors of J ftp decaying into pseudoscalar 
mesons could be regarded as the predictions of our model, as there are no experimental data 
at present. The introduction of light quark components in J/ip and rjc not only allows them 
to decay into the light mesons directly without intermediate gluons or virtual photon but is 
also helpful for us to understand the structures of charmonium states better. Considering 
that the mixing introduces a cc component into the light mesons, and such a cc component 
is intrinsic to the wave functions and exists over a time scale independent of any probe 
momentum, we could naturally interpret it as the intrinsic charm of these mesons. Our 
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Q'(GeV') 



FIG. 3: The Q^ behavior of the form factor Q'^ Fjji ^^^* (Q'^) compared with experimental data 42 
ii- 



result of the intrinsic charm content of the rj' meson fZ is also comparable with predictions 
from other models. 
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Appendix A: 

The SO (4) group elements can be written in terms of the SO (3) group generators [A^ 
and Bk): 

M = R+R_, (Al) 

i?+ = e-^^^-^S R_ = e-'^^+^^^ , {k = 1,2,3). (A2) 

The generators Ak and Bk obey the commuting relations of S0(3) generators [39|]: 

[Ai, Aj] = iSijkAk, [Bi, Bj] = iSijkBk, [Ai, Bj] = 0. (A3) 
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FIG. 4: Prediction of the Q^ behavior of the form factor Q'^Frj„^.yy*(Q'^), compared with the predic- 
ts). The dotted curve of the perturbative 



tions in the leading order of the perturbative approach 

approach indicates the Q^ region where QCD corrections may alter the predictions slightly. 



We see that the groups have the relation 5*0 (4) = SO (3) ® 5*0(3), and the generators A^ 
(as well as Bk) (k=l,2,3) could be seen as the angular momentum operators in each of the 
three directions. 
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FIG. 5: The Q^ behavior of the form factor Q'^F^^^^^* (Q^) compared with experimental data 46 
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From the matrix form of A^, we have 
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02 e, -01 
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An is the angular momentum component of the direction n. In fact, the matrix form of A, 
in Eq. (1A8I) can be diagonalized as 



a 2a 



A' 



( -\ o\ 
-i 



(A9) 



I 

V iy 

which is the expression of angular momentum operator An in its eigenstate representation, 
with the eigenvalues of An being (—1/2, —1/2, 1/2, 1/2). The two matrix forms are related 
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FIG. 6: The Q^ behavior of the form factor Q"^ F^p^y^^* (Q'^) compared with experimental data 48|. 



as A' = S^AnS, with the transformation matrix 



/ 0. 



S 



V2^ 



a 



02+id3a 










a 



8162— iOja 





a 



(AlO) 



In the eigenstate representation of angular momentum An, the matrix form of R+ is 

/ef \ 

et 

e"5^ 

\ e^ 



RL 



-iaAr, 



I 



and then we have 



/ cosf -^sinf 



R4. = SRj_S 



_92 



sm 



^sinf 

a 2 

02 



sm 



cos I 
^sinf 

a 2 

02 



-^sinf 

a 2 



COS ■ 



-^sinf \ 

a 2 

^sinf 

a 2 



^sinf 

a 2 



V^sin^-^sin^ ^ sin ^ cos ^ / 

\ a 2 a 2 a 2 2 / 



(All) 



(A12) 
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FIG. 7: Prediction of the Q^ behavior of the form factor Q'^ Fj/^^^^* (Q'^) 

which is the expression of R+ given in Eq. (j4]). Following the same procedure, we also obtain 
the expression of R- given in Eq. ([5]) . 



Appendix B: 

During the numerical calculation, we write M in a more compact form with eight real 
parameters (a, b, c, d, p, q, r, and s) under the constraints a^ + b"^ + c^ + (f = 1 and 
p^ + q^ + r^ + s^ = 1 J4l|: 



M 



f a -b -c -d\ ( 

b a —d c 

c d a —b 

\d —c b a J 

( 



p —q —r —s \ 
q p s —r 
r — s p q 
\s r —qp) 

ap—bq—cr—ds —aq—bp+cs—dr —ar—bs—cp+dq —as+br—cq—dp 

bp+aq—dT+cs —bq+ap+ds+cr —br+as—dp—cq —bs—ar—dq+cp 

cp+dq+ar—bs —cq+dp—as~br —cr+ds+ap+bq ~cs—dr+aq—bp 

\ dp~cq+br+as —dq—cp—bs+ar —dr—cs+bp—aq —ds+cr+bq+ap j 



(Bl) 



(B2) 
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FIG. 8: Prediction of the Q^ behavior of the form factor Q'^Fji^^^^* (Q^) 



These parameters are related to the six rotation angles as 



a 



Oo 



a 



a 



COS—, 0= — sm— , c= — sm— , a = — sm — 



P 



2' a 2 

COS-, g = -^sin- 



a 2 

On . P 
r = — - sm — . 

f3 2' 



a 

Oe . 13 
' = -13'''' 2- 



When referring to the mixing of specific types of mesons, the parameters (a, h, ... 
to (a^, 6^, ...) or (a^, hg, ...) correspondingly. 
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